Capture of particles of dust by convective flow 
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Interaction of particles of dust with vortex convective flows is under theoretical consideration. It 
is assumed that the volume fraction of solid phase is small, variations of density due to nonuniform 
distribution of particles and those caused by temperature nonisothermality of medium are compara- 
ble. Equations for the description of thermal buoyancy convection of a dusty medium are developed 
in the framework of the generalized Boussinesq approximation taking into account finite velocity 
of particle sedimentation. The capture of a cloud of dust particles by a vortex convective flow is 
considered, general criterion for the formation of such a cloud is obtained. The peculiarities of a 
steady state in the form of a dust cloud and backward influence of the solid phase on the carrier 
flow are studied in detail for a vertical layer heated from the sidewalls. It is shown that in the case, 
when this backward influence is essential, a hysteresis behavior is possible. The stability analysis of 
the steady state is performed. It turns out that there is a narrow range of governing parameters, in 
which such a steady state is stable. 
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I. INTRODUCTION 

Disperse systems, such as clouds of small particles sus- 
pended in liquid or gas are widespread in natural environ- 
ment as well as in various fields of human activity. Due 
to admixture, such systems demonstrate intriguing phys- 
ical effects, which cannot occur in homogeneous media. 
Understanding of the mechanisms governing the behav- 
ior of particles in fluid flows is of crucial importance both 
for fundamental studies and practical applications. 

Advection of particles in fluid flows is currently in- 
vestigated in several disciplines. In chemistry (see, for 
example, Refs. 00); this problem is essential in lower- 
ing fuel consumption during combustion of liquid fuels. 
In biology, recent developments have put forth a num- 
ber of fascinating problems concerning transport of mi- 
croorganisms in aqueous medium^*^ A typical example 
are algae, transported by ocean flows. These algae con- 
tribute essentially to the absorption of CO2 in oceans, 
which is closely related to the problem of climate change. 
In rapidly developing microfluidics, advection of particles 
has been the focus of close attention due to a need of find- 
ing effective mechanisms of particle mixing.® In fiuid me- 
chanics light-scattering particles of dust are widely used 
for visualization of flows. Besides, it is also important to 
know how particles influence the flow. Investigation of a 
spread of fine-dispersed impurities in the atmosphere and 
oceans is tightly related to the problem of environmental 
protection. These numerous applications have motivated 
the unfiagging interest in the fundamental problems of 
particle advection in fluids. 
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The behavior of a small single particle in isothermal 
fiuid fiows has been the subject of a great deal of study 
(see, for example, the monograph^ and the reviewsSiS). 
In particular, the equation of motion for a small single 
particle in a nonuniform unsteady fiow was discussed by 
Maxey and Riley>iS However, the fact that this equation 
takes into account the integral Basset force essentially 
complicates its solution. Therefore in most cases the 
problem has been solved in terms of asymptotic theo- 
retical models. Stommelii studied sedimentation of par- 
ticles in a cellular fiuid fiow in terms of weak inertia ap- 
proximation. To our knowledge this was the first paper 
that developed a simple description of capture of small 
solid particles by a fiuid fiow. It was shown, that when 
the ascending fiow is rather intensive, some particles are 
involved in a vortex motion and remain suspended. Re- 
cently, capture of particles has been observed in a cel- 
lular convective fiow.-i^ Qualitatively similar behavior is 
shown by granular media in the ascending air fiow un- 
der gravity: at a certain velocity of air fiow the granular 
matter is liquidizedi^^ The problem of sedimentation of 
small heavy particles has been generalized for the case 
of a random fiow (see Ref. ^3 and references therein). 
However, a backward effect of the particles on a fluid 
flow was not considered. A number of theoretical and 
experimental studies are devoted to the problem of par- 
ticle accumulation in fluid flows - the effect, which is 
essential due to a difference in the inertial properties of 
the phases when gravity is insignificant ii^*iSiiLi& The sit- 
uation, when both gravity and inertia contribute to the 
particle dynamics, has been studied in papers li^iSSiSi The 
authors of papersSSiSi^^iSiSSiSi have investigated the in- 
fiuence of the Basset history force on the particle motion. 

Another problem that has been intensively studied 
for many years is the stability of laminar isothermal 
fiows carrying small solid particles. First papers22i22i2£ 
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were devoted to the stability of the particulate Poiseuille 
flow. This problem has been generalized in the series of 
papers i?^i?^ The stability of the Couette flow was consid- 
ered in Ref. 33. A general result of these investigations 
is that stability of a dusty flow qualitatively depends on 
the size of particles: relatively small particles destabilize 
the flow, whereas the particles of a larger size make the 
flow more stable. The relaxation time of small particles 
is short, therefore such particles move actually with the 
velocity of the fluid. The presence of particles results in 
renormalization of the media density. A larger efficient 
density is equivalent to a higher fluid velocity, which leads 
to destabilization. On the contrary, the relaxation time 
of large particles is long. Such particles are too inert to 
respond to rapid velocity variations of the carrying fluid 
and hence, large particles damp the fluid velocity pertur- 
bations. 

The collective behavior of particles in laminar non- 
isothermal flows is much less studied. Moreover, the 
problem of interaction of the particles with convective 
flows is not completely understood yet. The exceptions 
are the works, concerning the influence of settling parti- 
cles on the stability of convective flowsi^^i^SiSfe 

It is known, that particles of dust are widely used for 
visualization of flows in the experiments on thermal buoy- 
ancy convection. It is conventionally assumed that, if the 
mass concentration of particles is small, their backward 
influence on the convective flow is negligible. However, 
this is not always true. The point is that buoyancy con- 
vective flow itself is caused by small variations of density 
related to nonisothermality of a medium. If density vari- 
ations due to nonuniformity of the particle distribution 
are of the same order of magnitude as those due to non- 
isothermality, one should no longer ignore the influence 
of particles on the flow. For small particle concentration 
this effect is intuitively believed to be pronounced if the 
density of particles and the density of the carrier fluid 
differ enough. A typical example is solid particles sus- 
pended in a gaseous medium. However, we argue below 
that the effect is essential in more general case, when the 
densities of phases are different, but of the same order 
(for example, in liquids laden with small solid particles). 
Under gravity such particles tend to settle on the bottom 
of a container. Two mechanisms hinder particles from 
settling down: the Brownian motion of particles leading 
to their diffusion in a carrier fluid and the entrapment of 
particles by vortex fluid flows. The first mechanism op- 
erates also in a quiescent medium and leads to the Boltz- 
mann distribution of particles over a height. However, 
for particle sizes commonly used for flow visualization, 
an effective diffusion coefficient is so small that in the 
equilibrium state the particles are practically absent in 
most part of the volume. The second mechanism, which 
operates only if the maximal velocity of ascending flows 
exceeds the sedimentation velocity, is much more effi- 
cient. Investigation of the capture effect and a backward 
influence of the particles on a flow is important not only 
for visualization problems but also for more fundamental 



aspects of environmental problems. 

The mechanism of particle capture discussed above is 
related to nonuniformity of the particle distribution over 
the volume. Even if the initial distribution of particles 
is uniform, this nonuniformity arises by itself due to the 
particle sedimentation in the regions, where the flow is 
either descending or has small vertical velocities. If the 
vertical dimension of the container is large, the particle 
distribution in the most part of the container is nearly 
uniform for a long time. However, even in this case the 
particles can influence the flow, because their velocity dif- 
fers from the velocity of the fluid due to sedimentation. 
Thus, the particles can transfer the perturbations of tem- 
perature and vorticity and affect stability of a convective 
flowM 

In the present paper, we study the capture of dust par- 
ticles by convective flows. Particularly, a backward effect 
of particles on the fluid flow is investigated in details. In 
Sec.^ an appropriate theoretical model is developed as a 
generalization of the Boussinesq approximation to a con- 
vective flow in a dusty medium. In Sec. lIIII the obtained 
equations are used to make a general interpretation of the 
problem and to prove that the entrapped particles form 
a cloud of dust. Section El is devoted to consideration 
of a dusty medium in an infinite vertical layer. The exis- 
tence of a steady state solution with particle distribution 
in the form a dust cloud is proved, and the stability of 
this state is investigated. 



II. THEORETICAL MODEL OF THERMAL 
BUOYANCY CONVECTION IN A DUSTY 
MEDIUM 

A. Governing equations and basic assumptions 

Consider the behavior of small solid particles, sus- 
pended in a nonisothermal fluid (liquid or gas) under 
gravity. On one hand, all particles are assumed to 
be monodisperse spheres of a radius rp, which is large 
enough to neglect particle diffusion. On the other hand, 
this size is supposed to be much smaller than the char- 
acteristic length scale L of the flow. On a scale much 
greater than the particles are regarded as a continu- 
ous medium with the volume fraction ip = A/Snrpn (the 
volume fraction of the fluid phase fi = 1 — tp), where 
n is the the number of particles per unit volume of the 
medium. Since the actual volume concentration n is pro- 
portional to if and differs from the latter by a constant 
factor, in the following for the sake of simplicity if will 
be called concentration. 

We suppose that the volume concentration of parti- 
cles is small, so we can neglect interparticle interactions 
and interactions between the particles and walls of a con- 
tainer. It is also assumed that the carrier phase is in- 
compressible and the solid particles can neither deform 
nor combine into agglomerates. Moreover, the density 
of a solid phase is considered to be constant. The latter 
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assumption implies that we neglect thermal expansion of 
the solid phase, which is perfectly valid in the dilute limit. 
After averaging over space the equations for mass, mo- 
mentum and energy balance of both phases are written 
as followstiSiSi 



dt 



+ div (/ipu) = 0, 



— + div ((^Up) = 0, 



(la) 
(lb) 



^p— = -Vp + V • {^le) - fipG&e^ - tpF, (2a) 

= -(5Gae, + F, (2b) 

DT 1 3 
MP— - — div(MVT) + — — ^(Tj,-T), (3a) 

Pr rf' 



Dt Pr 



p 



at Pr r£ 



(3b) 



where u and Up, T and Tp, are the velocities and the 
temperatures of phases, p and p are the pressure and 
the density of fluid, respectively (hereafter, the subscript 
"p" stands for the particle phase), is the unit vector of 
the axis z, directed against gravity; the shear rate tensor 
Eij = S/iUj + \/jUi. The dimensionless variables have 
been introduced using the following units: the reference 
density of fluid po for the densities of phases, the reference 
temperature difference 9 for the temperatures, L for the 
coordinates, L^/v for the time, ly/L for the velocities, 
Pqv'^ / for the pressure, where ry is the viscosity of fluid 
and V = rj/ pq. 

Equations Q-® involve the following dimensionless 
parameters: the Galilei number Ga = gL^/v'^, the 
Prandtl number Pr = rfc/K, the ratio of densities of 
phases 5 = ppj p^^ the ratio of the specific heats of phases 
at constant pressure B = Cp/c and the ratio of character- 
istic particle size to the flow length scale r-p/L, where k, 
is the thermal conductivity of fluid. We make a note of 
the distinction between the two Lagrangian derivatives 
DjDt = a/9t-|-u-V and djdt = d/dt + Up-V, which are 
used to denote the time derivative associated with the 
motion of the fluid element and the element of a solid 
phase, respectively. 

In the interphase interaction term, F has the meaning 
of the force exerted by an unsteady nonuniform fluid on 
the solid particle. This force can be written asiSiS 



= P 



Du 

Ijt 



9 L2 

Gae.)-^-W 



9 L 
2^ 



dW{T) dr 
dr y/T^ 
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where the relative velocity of phases W = Up — u is in- 
troduced. The first two terms in I^J , correspond to con- 
tributions to the force exerted on the particle by undis- 
turbed fluid flow due to gravity and pressure gradient. 



In the approximation of undisturbed flow these terms 
coincide with the Archemedian force. Generally, the par- 
ticle during its motion, disturbs the flow, the effect being 
taken into account by the next three terms in Q): the 
Stokes viscous drag, the Basset history force, caused by 
unsteadiness of a viscous boundary layer around a parti- 
cle, and the added mass force, allowing for the inertia of 
the surrounding fluid. 

We do not consider the effects that can be initiated by 
rotation of particles. Particularly, the Magnus force is 
not taken into account in Q). This assumption holds only 
in the case when the particle response time 2r'^5/{%v) 
is much less than the characteristic hydrodynamic time 
scale /v. Therefore the Stokes number, corresponding 
to a ratio of these time scales, is assumed to be small: 

9 

However, we assume the parameter 6 to be finite and do 
not restrict our theory to the limiting case of dusty media 
(<5> l),cf. RefsT^T?. Further, we do not explicitly take 
into account Einstein's correction term to the viscosity 
of fluid due to solid admixture, which is inessential in the 
dilute limit. 

The equation of state used here is typical for ther- 
mal buoyancy convection. We consider weak convection, 
when the variations of density due to nonisothermality 
of a fluid are small. Assuming that the fluid density is 
a function of only temperature p = p{T), we expand it 
into a Taylor series near its reference value po at the 
temperature value Tq and restrict ourselves to a linear 
approximation 

p=l- per, 139 < 1, 

where temperature T is measured from the value Tq, /? 
is the thermal expansion coefficient. 



B. Single-fluid approximation 

Let us consider a fluid laden with small solid particles 
and prove that the initial two-phase theoretical model 
(^-Ql can be simplified. We assume that the mass con- 
centration of particles dip is comparable with the relative 
variations of the fluid density due to nonisothermality. 
Bearing in mind that the particle concentration ip, the 
nonisothermality P9, and the relative size of particles 
Tp/L are small we retain only the leading terms in the 
momentum equations l(2Jl. As a result we conclude that 
in the basic state the pressure distribution is hydrostatic, 
and obtain a relation between the velocities of phases: 



Vp-I-Gae^ = 0, 



(5) 



u - Se^, 



2^(^_l)Ga, (6) 



where S is the dimensionless sedimentation velocity. 
Equation (0 indicates that in the accepted approxima- 
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tion the particle velocity equals the velocity of the fluid 
plus the constant sedimentation velocity. 

Taking into account we obtain in the next order 
the momentum equation of the fluid 

-^ + u-Vu^ -Vp'+V^u+i/SeT - {S - l)ip]Ga.e,, (7) 

where p' is the convective addition to the hydrostatic 
pressure. Following the idea of the Boussinesq approx- 
imation we can state that at large values of the Galilei 
number and small density inhomogeneities the density 
of the medium will be constant everywhere except for 
the buoyancy force, where these small density variations 
are multiplied by Ga. So, we introduce the thermal 
Grashof number Gr = (39Ga and its concentration analog 
Gc = ((5 — l)(poGa, where ipo is the characteristic value 
of the particle concentration. 

From the energy balance equation Ij3b|l we obtain in the 
leading order: Tp = T. We restrict our consideration to 
this approximation. This implies that for small particles 
the time to equilibrate the phase temperature is much less 
than the viscous hydrodynamic time scale. The energy 
equation for the fluid phase is then obtained from Eqs. 
®. 

Taking into account the mass balance equations for 
fluid and particles ^ we arrive at a complete set of 
equations describing thermal buoyancy convection in a 
nonisothermal fluid laden with particles: 



du 

'dt 
dT 

'dt 

dip 



H u Vu 

-u - VT 

Up • \/ip 
div u 



-Vp + V^u + (Gr T - Gc(p) e^, (8) 



1 

ft 

0, 



= 0, u„ = u — Sgz 



(9) 

(10) 
(11) 



where the prime for p is omitted, and the concentration 
is normalized by ipo- 

The set of equations (|5|)- (tTT|l is reminiscent of equa- 
tions for thermosolutal convection. However, there are 
two essential differences: first, we have neglected the dif- 
fusion of particles, second, we have taken into account 
the finite sedimentation velocity. The developed model is 
quite general and can be applied to dusty media, aerosols, 
liquids laden with small solid particles, and biological 
species in aqueous media. 



III. CAPTURE OF DUST PARTICLES. 
GENERAL CONSIDERATION 

Consider the behavior of a dusty medium in a closed 
cavity of an arbitrary shape heated from the sidewalls. It 
is assumed that at the initial time the dust is uniformly 
distributed all over the volume and at a later time no new 
particles enter the cavity but the existing particles may 
settle to the bottom. We also assume that the particles 



that have settled on the bottom are not carried back into 
the flow, but stick to the bottom. The question arises as 
to whether all particles should eventually settle on the 
bottom or there may occur a situation when some parti- 
cles will be entrapped by the flow and stay in a suspended 
state for inflnitely long time. 

Obviously, in the absence of heating (Gr = 0) and 
any initial perturbations, all particles will eventually sink 
down on the bottom for a time of order 1/S, which corre- 
sponds to the following solution of the equations (|S|)- (|11|I : 



u 

T 



0, z> Z{t) 

1, z < Z{t) 



where Z{t) — const — St. So, initial distribution of par- 
ticles ensures sedimentation without distortion of the in- 
terface between a pure fluid and suspension. 

In order to simplify further argumentation, we will 
measure z from the lowest point of the cavity. It can be 
understood through insight that inducing of initial per- 
turbations of concentration and velocity does not change 
the situation: finally all perturbations will be damped 
and all the particles will settle on the bottom. To demon- 
strate this let us consider the evolution of the total energy 
of the system E: 



E 



Gc tfz I (Kl 



(12) 



Using (|5|l- Hll|l we obtain for the rate of the energy 
change 



= — / (V X ufdD, + i) ipze^ -n dF, 
at Jn Jp 



(13) 



where F is the surface bounding the volume fJ, n is the 
external vector normal to F. To obtain Eq. (|13|l the no 
slip conditions for fluid velocity are assumed: u\p = 0. 
Since at all time moments t > the value (p differs from 
zero only on such areas of the boundary surface where 
ez • n < (i.e. at the bottom of a cavity), the right-hand 
side of H13() is always nonpositive. Moreover it can be 
equal to zero only in the case of a quiescent fluid in the 
absence of the particles at the walls. This proves that the 
energy E can only decrease. On the other hand, as can 
be seen from (|12|l . E >0, i.e. it is bounded from below. 
Hence, the flnal state is that, at which the right-hand 
side in (|13|l turns to zero, i.e. when the fluid becomes 
completely free of the particles and all fluid flows cease. 

The situation radically changes if there is a permanent 
external source that makes fluid move. In the present 
work, the source of this kind is side heating of the wall, 
generating a vortex convective flow in a cavity. Naturally, 
in this case, the areas near the walls will finally be free 
of particles. Indeed, as it follows from Eq. H1U|) 



dtJn 



Up ■ \7ip dfl ~ S (p ipe^ ■ n dF < 0. 



(14) 
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FIG. 1: Streamlines of the particle flow in the absence (a) 
and in the case (b) of particle capture. 

The total amount of the dust in the cavity decreases un- 
til <f turns to zero everywhere near the boundaries. It 
does not mean, however, that the bulk of the fluid will 
be completely free of particles. If in some part of the 
cavity the velocity of the ascending fluid Ua exceeds the 
sedimentation velocity S, we can expect the capture of 
the particles. The condition 

Ua > S (15) 

is the necessary condition for capture but it is not a suffi- 
cient one. Indeed, we can easily imagine the flow pattern 
with the streamlines of particles shown in Fig.^a): de- 
spite the existence of regions with ascending motion of 
particles, all streamlines end at the lower boundary. 

On the other hand, we can specify a sufhcicnt condi- 
tion for capture, which is probably not necessary. Let 
us assume that there is the line, at which the horizon- 
tal component of the the fluid velocity vanishes. If on 
this line there is a point where Ua = S then such a point 
is a singular point for the vector field Up. Furthermore, 
owing to 1)11(1 we have 

div Up — 0, 

and therefore the vector field Up is free of sources or sinks. 
In this case only two types of structurally stable fixed 
points are possible: center or saddle. In the case of the 
center point, there are closed orbits, i.e. some particles 
are entrapped by a vortex and stay inside for infinitely 
long time. The existence of the saddle point together 
with the boundary conditions for Up necessary implies 
the existence of the closed separatrix loop and the center 
point, as it is schematically presented in Fig. ^b). The 
separatrix loop forms the boundary of the captured cloud 
of particles. 

To demonstrate how formation of a cloud of dust oc- 
curs, we have numerically integrated Eqs. (|5|)- (|ll|l for a 
partial case of the square cavity, heated from the right 
side. The quiescent state with uniform distribution of 
particles was chosen as the initial one. The values of gov- 
erning parameters correspond to the single-vortex flow, 
satisfying the condition H15|l . This complementary exam- 
ple is in agreement with the results of general considera- 
tion. During evolution (see Fig. [21), a part of particles is 



gradually leaving the flow, whereas other particles stay 
suspended. As a result of this transient process, the sys- 
tem evolves to the steady state with a distribution of 
particles in the form of a cloud. 

The theory above describes the case of no particle in- 
ertia. In general, the trajectory of a captured particle 
is no longer a closed orbit. Due to inertia the center 
fixed point [see Fig. ^b)] transforms to a focus, which 
is unstable for heavy particles.^ In the vicinity of the fo- 
cus, captured particles of dust are repelled by this fixed 
point and move along spirals. Nonetheless, for small par- 
ticles this effect is weak. The addition to the velocity 
of particle due to purely inertial drift u™ is giveniSiii 
by the term StDw./ Dt. Since the intensity of thermal 
convection is governed by the thermal Grashof number, 
we can estimate u^" ~ St Gr. For fine particles of dust 
with Tp ~ 10~^ cm suspended in air on laboratory length 
scales L ^ 1 cm at moderate intensity of thermal con- 
vection we obtain the inertial time scale Tin ^ 10^ s 
(St - 10-5, Gr - 103). 

Recall, that we have also neglected the diffusion of the 
particles, which means that our consideration concerns 
the time scales less than the diffusion time scale td- How- 
ever, this restriction is practically of no significance. In- 
deed, according to the Einstein formula, the diffusion co- 
efficient of the spherical particles D = kT / {Gnrirp) (here 
T is the absolute temperature, k is the Boltzmann con- 
stant). As before, for the fine particles of dust in air at 
room temperature we obtain D — 10"^ cm^/s, which for 
the laboratory length scales corresponds to the diffusion 
time scale td ~ 10^ s (rp — lO""' cm, L — 1 cm). 

Thus, despite the fact that existence of a state with a 
cloud of dust is in general not everlasting due to particle 
inertia and particle diffusion, it exists long enough to be 
important. 




FIG. 2: Evolution of particle concentration (fi to the steady 
state at Gr = 10 000, S = 5, Gc = 0. The times for the states 
of (fi are: (a) t = 0.02, (b) t = 0.04, (c) t = 0.07, (d) t = 0.09, 
(e) t = 0.14, (f) t — 0.40. The regions of pure liquid and 
particle suspension are shown in white and grey respectively. 
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IV. CLOUD OF DUST IN A VERTICAL LAYER 

A. Basic steady state 

A rather complete understanding of the mechanisms 
responsible for formation of a cloud of dust and its back- 
ward influence on the flow hydrodynamics can be ob- 
tained from a model problem. Let a dusty medium fill 
the infinite vertical layer — f < a; < f , the boundaries of 
which are kept at constant, but different temperatures: 
T = — lata; = — 1 and T = 1 at a; = 1 . We assume that 
initial distribution of particles is uniform and look for a 
steady state solution, in which the velocity has only ver- 
tical component uq and all fields except for the pressure 
are independent of the vertical coordinate z. Then, the 
set of equations takes the form: 

Uq + Gr To — Gc(^o = c, (16) 
W = 0, (17) 

where the prime stands for differentiation with respect 
to X, the subscript "0" is used to indicate the steady 
state solution, c = dp/dz — const. Equations l(T7|) 
should be supplemented with the boundary conditions 
for the fluid velocity and the temperature: 

X = ±1 : uo = 0, Tq = ±1. (18) 

In the case of the infinite layer we should also specify 
integral conditions. We prescribe no-fiux for fluid and 
particles, which mean that the flow is closed at the infin- 
ity and no new particles enter the system: 

1 1 
JuQtlx — O, J (poUpodx = 0, (19) 
-1 -1 

where Wpo = — S. 

According to 117|) . (|18|l . the temperature distribution 
does not depend on the fluid and particle motions and 
can be obtained directly 

To = X. (20) 

Let us analyze velocity profiles of the fluid and solid 
phase. We start the discussion with the simplest case 
Gc — 0, when the particles do not influence the fluid 
motion. In this case, the velocity profile of the fiuid is 
exactly the same as in the well-known problem on a con- 
vective flow of pure fluid in a vertical layer heated from 
the sidewallsi^ 

uo = ^x{l~x^). (21) 

The particles occupy the region of the layer between 
Xi and X2- Within this interval ipQ ~ 1, whereas in the 
near- wall regions ipo — 0. 
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FIG. 3: The width of a dust cloud d versus Gro at Gc = 0. 

The point X2 is the point at which the velocity of par- 
ticle sedimentation coincides with the velocity of the as- 
cending fluid flow 

X — X2 ■ Uq — S. (22) 

According to H21|) , it is determined by the largest root of 
the cubic equation 

Gra;2(l-2:^) =6S. (23) 

The point xi is determined by the zero-flux condition 
(|19|l for particles. It is convenient to rewrite the equation 
for xi in terms of semi- width of the cloud d = {x2 — xi)/2 
(hereafter, just width). Thus, we obtain for d 

d = X2-\[^^. (24) 

As evident from Eq. 123|) and the second condition in 
(|19|l . the parameter S is not independent; the width of the 
cloud is determined by the ratio Gro — Gr/S, moreover, 
the cloud exists only if Grp > 9-\/3. With the growth of 
Gro the width of the cloud monotonically increases and 
tends to 1 as Gro — > oo (see Fig. I^)). 

At nonzero values of Gc the particles exercise an in- 
fluence on the fluid flow, and the dependence of the 
characteristics of the cloud on the parameters becomes 
more complex. In this case the basic state is gov- 
erned by Gro and a renormalized concentration parame- 
ter Gco = Gc/S. In the presence of a particle cloud the 
flow region splits into three zones: 1) — 1 < a; < xi, 2) 
xi < a; < X2, 3) a;2 < a; < 1- In each of these zones 
the velocity profile is described by third-degree polyno- 
mials, which can be written down in the following form, 
respectively: 

Uq^^ ~ uo + c (l — x^) + ci (1 + a;) , 

Uq^^ — uo + - Gco + c (l — a;^) -I- C21 a: -I- C22,(25) 

Uq^^ = Uq + c (l — a;^) -I- C3 (1 — a;) , 



7 




where uq is given by (|21|l normalized by S; the unknown 
coefficients of the polynomials ci, C21, C22, C3, the pa- 
rameter c and the coordinates of the points xi and X2 
are defined by the boundary conditions (|18|l . H19|l . the 
condition (|22l) and the continuity conditions for the ve- 
locity and tangential stress at the cloud borders. The 
system of the obtained algebraic equations is partly sim- 
plified: the constants c, ci, C21, C22, C3 can be excluded by 
expressing them in terms of the width of the dust cloud 
d and the coordinate of its center Xc = {x2 + Xi)/2: 

c = dGco [3{l~xl)-d^] , 
ci = -dGco (1 - a;c) , C3 = -d Gcq (1 + Xc) , 

C21 = -Geo Xc (1 - d), C22 = ^ Geo (xl + d^ - 2d) . 

The remaining unknowns d and defined by the 

couple of nonlinear algebraic equations, which are solved 
numerically 

12 = 3dGco{xc + d-l)fi{xc,d) 

- 2 Gro {xc + d) Uxc + df + ll , (26) 
12 ^ 3dGcof2ixc,d)~2GToXc{xl + d^ -l),{27) 

fi(xc,d) = 3xl + 3{d+l)xl + {d'^ + l)xc 

+ d^ + d^ - 3d + I, 
f2{xc,d) = {3xl + d^ + 2d-3){3xl + d^ -2d+l) . 

As in the case discussed above, we assume that Lp^ — Q 
in the zones 1 and 3 and Lp^ = \ in the zone 2. 

The cloud is absent if the velocity of fluid is lower than 
the velocity of sedimentation all over the volume. As in 
the case of Gcq — 0, this means that the cloud can exist 
only at Gro > 9-\/3. 

Let us discuss the dependence of the width of the par- 
ticle cloud on the concentration parameter Gco- At a 
relatively low Gro the width of the cloud d monotoni- 
cally decreases with the growth of Gcq (see Fig. 0] line 



1) . This is related to the fact that the major part of the 
cloud is located in the ascending flow, hence, the larger 
Gco, the higher effective density of the medium in this 
place. The higher density leads to a decrease of the buoy- 
ancy force and therefore to lowering of the flow velocity. 
The latter, in its turn, results in a decrease of the width 
of the dust cloud, in which the particles can remain in 
the suspended state. 

At higher values of Gro the dependence of d on Gco be- 
comes more complex. Generally, in this case the growth 
of Gc also results in the decrease of the cloud width, but 
now in some range of Gco the dependence (i(Gco) is no 
longer unique (see Fig. 01 line 2). At low Gco the growth 
of Gco only slightly influences the width of the cloud, 
but starting from some threshold value of Gco, the drag 
force of the wide cloud becomes so strong that the flow 
is unable to keep it further. As a result, the width of 
a cloud decreases by a jump. If now, starting from a 
large value of Gco, we decrease it, first, the cloud width 
varies insignificantly, but then at some critical value of 
concentration parameter the width of a cloud makes a 
jump to nearly the highest possible value for a given Gro. 
Thus, the decrease of the cloud width with the change of 
particle concentration demonstrates a well pronounced 
hysteresis. 

In Fig.|Slthe profiles of the fluid velocity at Gro — 500, 
Gco = 600 are plotted for upper (Hue 1) and lower (Hue 

2) branches of possible steady states. It is seen, that for 
the upper branch the velocity profile only slightly differs 
from the usual cubic profile (|21|l , but for the lower branch 
the flow in the domain occupied by the cloud is strongly 
suppressed, the total intensity of the flow is also low. 

In Fig. we plot a diagram defining a range of hys- 
teresis existence on the plane (Gco, Gro). The solid lines 
divide the plane into two parts: larger and smaller ones. 
Any point of the larger part corresponds to only one so- 
lution (one root of a cloud width d at fixed Gco, Gro), 
whereas in the smaller one, the hysteresis zone, there are 
three solutions (see Fig. The solid lines, at which 
there exist two solutions, correspond to folds; the point 
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FIG. 6: Hysteresis range on the parameter plane (Geo, Gro). 



where these hnes intersect is a cusp with one solution. In- 
side the hysteresis zone, one of three solutions, namely, 
corresponding to the middle branch of a hysteresis curve 
(i(Gco) in Fig. ^ is structurally unstable. 



B. 3.2. Linear stability analysis 

1. Formulation of stability problem 

Let us investigate the stability of the basic steady state. 
We restrict our consideration to a fixed value of the 
Prandtl number Pr = 1, which corresponds to a typi- 
cal case of particles, suspended in gaseous medium. 

In order to formulate the stability problem, the govern- 
ing equations lP^- Hll(l should be supplemented by neces- 
sary conditions at the rigid boundaries of the layer and 
continuity conditions at the borders of the dust cloud. 
At the rigid walls of the layer we impose the no slip con- 
dition for fluid velocity, and the condition of constant 
temperature: 



x = ±l: u = 0, T = ±l. 



(28) 



The borders of the dust cloud, which are the surfaces of 
discontinuity for particle concentration, must satisfy the 
conditions of stress balance, the continuity of velocity, 
temperature, and energy flux. These conditions are the 
direct consequence of conservation laws for momentum, 
energy, and mass. Assuming that the borders of a dust 
cloud are deformable, the conditions at the interfaces are 



X = Xi^2 + Cl,2 : + Km] = 0, [Snr] = 

'dT 



[u] = 0, [T] = 0, 



dn 



0, (29) 



where Ci = C,i{y, z, t), C,2 — C2(y, t) are the deviations of 
concentration surfaces from the flat undisturbed shapes, 
and square brackets are used to denote the jump of a 
function [/] = /2 — /i,3- Here the subscripts 1, 2, 3 
correspond to three zones of the flow: 1) — 1 < a; < 



xi + Cij 'i) xi + Ci < X < X2 + C2, 3) a;2 -I- C2 < a; < 1. 
The unit normal n and tangential Ti, T2 vectors to the 
surface are given by the relationships 



n = 



n = 



'l + (VCl,2)' 

r2 = n X ri . 



; + Bx ■ VCl,2 

1 + (e. • VCi,2)' 



Here we introduce the orts Bx — (1,0,0), e^, = (0,1,0) 
of the axes x and y, respectively. Both borders of the 
dust cloud are impervious to the particles, the velocity 
of any point of such border coincides with the Lagrangian 
velocity of a particle at this point. Hence, for surface de- 
viations the following kinematic conditions are specified: 



1,2 



dt 



u„ 



VC: 



1.2 



(30) 



Let the basic steady state be disturbed by introducing 
small perturbations of the velocity v, pressure g, tem- 
perature "d and concentration 0. Then we substitute the 
disturbed fields Uq-I-v, po + q, Tq + i), (po + 4' into Eqs. ©- 
(|ll|l . Neglecting squared and higher order terms with re- 
spect to perturbations, we obtain the equations for per- 
turbations: 



— uo • Vv - 

ot 



— +uo-V'd 
ot 



V • Vuq 



Se, 



VTo 



-Gc(t> 



= 0, 



V^v + Grz^e^ 
e., (31) 
, divv = 0,(32) 

(33) 



where Uq = (0, 0, uo(a;)). 

The boundary conditions H28|l take the form for per- 
turbations: 



±1 : v = 0, d^O. 



(34) 



Let us also assume the smallness of Ci, C2- Then, the 
equations H31|l - H33|) . boundary conditions H34|l and con- 
ditions at deformable interfaces H29|) reduced to those for 
undisturbed interfaces admit a transformation, which is 
analogous to the Squire transformation"^^ and discussed 
below. Such a transformation allows us to reduce the full 
three-dimensional problem to a problem in two dimen- 
sions. We denote the x and z components of the velocity 
V by M and w, respectively and analyze the stability of the 
basic steady state with respect to the two-dimensional 
perturbations in the form of transversal rolls. Thus, the 
solution can be written as normal modes 



/ u{x,z,t) \ 
w{x, z, t) 
q{x,z,t) 
d{x,z,t) 
<p{x,z,t) 

V Cl,2(^,t) / 



/ n{x) \ 
■w{x) 
q{x) 
{){x) 

\ Cl,2 J 



^Xt—ikz 



(35) 
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where A = + iXi is the complex growth rate, k is the 
wave number, characterizing periodicity of perturbations 
along the z-axis. 

Let us substitute the ansatz (|35|l into Eqs. (|31|) - (|33|l . 
It follows from H33() that 

{A-i/c(wo-S)}0 = O, 

which must be satisfied at every point of a layer occupied 
by particles. Since uq is a function of x we assume that 
{X — ik{uo — S)} ^ 0, (j) — 0. As a result, in each of 
the three zones the problem is defined by the same set of 
amplitude equations: 

Xii ~ ikuou = —q' + u" — k'^ii, (36) 
Aw — ikuoii + u'qU = ikq + w" — k^w + Gr?9, (37) 
u-ikw = 0, (38) 

Ai? - iku^d + u = —{&'- k^d] . (39) 
Pr V / 

After linearization, the boundary conditions (|34|l and 
the conditions at the pure liquid-suspension inter- 
face, i.e. at the borders of a dust cloud, written for the 
amplitudes of perturbations, take the form: 

X = ±1 : li = 0, ui = 0, = 0, (40) 
X = xi,2 ■■ [q] = 0, [it] = 0, [w] = 0, 

K]--GcCi,2, [^]=0, 0,(41) 

where C^i and C,2 are obtained from Eqs. (|30|l together 
with H41|) and are determined by the relations 

X = xi : {A — ik{uQ — S)} Ci ~ u] x ~ X2 '■ XC2 = u. 

(42) 

The conditions H40|) - (|42|l are defined taking into account 
the following properties of the undisturbed velocity pro- 
file (EHJ: [uo] = 0, K] = 0, [<] = Gc, and uq{x2) = S. 

The boundary value problem is a spectral 

amplitude problem. The conditions of nontrivial solution 
define perturbation growth rate as a function of the pa- 
rameters Gr, Pr, Gc, S and k. The solution with A^ = 
determines the neutral behavior and separates the re- 
gions of unstable modes with A^ > from those of stable 
modes with A^ < 0. 

In the particular case of Gc = 0, the particles have no 
influence on the fluid motion, and the boundary value 
problem H3()|l - H42l) reduces to the stability analysis^ of a 
flow with the odd velocity profile (|21|l . For Pr = 1, the in- 
stability appears above the critical value Gimin = 496.3, 
which is reached at kmin = 1.404, and corresponds to 
the nionotonic (i.e. to the solution with zero imaginary 
part of the growth rate A^ = 0) "hydrodynamic" pertur- 
bations. 

In the case of Gc ^ 0, the problem admits analyt- 
ical solution in the limit of short wavelength behavior 
for a small, but finite width of the dust cloud. In the 
general case, the boundary value problem 1)361) - (|42l) was 
treated numerically by the standard shooting and differ- 
ential sweep methods, which gave very close agreement 
of the results. 



2. Thin dust cloud, short wavelength limit 

In the case of a thin cloud, the parameter d is small, 
and the problem can be simplified. The nonlinear al- 
gebraic equations lf77|) . defining the width of the 
cloud d and the coordinate of its center be solved 
explicitly with the accuracy O (d^) : 

In the limit of vanishing cloud width, when d — > 0, we 
obtain the lower threshold for the existence of a dust 
cloud Gr = 9\/3S, which is in agreement with the the 
value obtained above for an arbitrary value of d. It is 
clearly seen from (I25f) that in this case the velocity profile 
coincides with (|21|l . A thin cloud is located in the vicinity 
of the point at which the fluid velocity is maximal. 

The short wavelength limit means that the wave num- 
ber is large fc ^ 00. Let us investigate the case of small, 
but finite values of d. Formally, this case corresponds to 
the limiting transition, when ^ 0, fc — )■ cx3, but their 
product kd is finite. 

This particular case can be examined in the context of 
the auxiliary problem, in which the layers of pure fluid 
are semi-infinite and the dust cloud is in between these 
layers and has the finite width 2d. It is convenient to 
treat the problem using the multi-scales technique4& Let 
us introduce the "fast" coordinate ^ = (x — Xc)/s, where 
£ — 1/fc is the small parameter. In terms of the fast 
coordinate the dust cloud occupies the finite region — A < 
^ < A with the center at the point ^ = 0. Here A — kd 
is a semi- width of the cloud, measured in the units of the 
fast coordinate. The layers of pure fluid occupy the areas 
— 00 < ^ < — A and A < ^ < cxd, respectively. 

The boundary value problem H36l) - (|42|l is rewritten in 
terms of the fast coordinate ^. As a result we obtain the 
equations: 

e^Xu — ieuQU = —eq' + u" — u, (44) 
e'^Xw — ieuQii -\- s^u'qU — ieq + w" — w 

-fe^Gri?, (45) 
u'-iw = 0, (46) 

e^Xd-ieuo^ + e^ii = ^ i^" ~ , (47) 
Pr V / 

which are the same for all three zones: 1)— cx)<^<— A, 
2) -A < ^ < A, 3) A < ^ < 00. We require that the 
functions ti, w, q, 1) must be bounded at ^ = ±00. The 
conditions at the borders of the cloud arc as follows: 

C ±A : [q]^ 0, [u] ^ 0, [w] = 0, 

H = -£GcCi,2, m^Q, m= 0,(48) 

C = -A: {£A-i(Mo-S)}Ci (49) 

^ = A : AC2 - ^i, (50) 

where Ci and (^2 refer to the points ^ = — A and ^ = A, 
respectively. 
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For each zone the solution is sought as a power series 
in £. Taking into account l|43|l . we rewrite the velocity 
profile (|25|l with the accuracy O (e'^) : 



,(1) 



,(2) 



,(3) 



S-.4S(ii# + Af^-Gc)(e + A) 



lx/3" 

S-eU -^li^^ - A ( ^ - Gc ) (e - A) 



V3 



Gr 



The boundary value problem 144() - (|50|l is the inner one 
with respect to the initial problem H3t)|l - H42|) . which in its 
turn should be referred to as an outer one. Generally, 
the solutions of inner problem is matched to the outer 
one. However, this should not necessary be done in the 
present case since the solution of the inner problem H44|) - 
(|5n)l proves to be vanishing at ^ = ±oo. 

The solution of the boundary value problem l|^ - l(Sn|) 
is quite cumbersome and is not adduced here. Actually, 
the expressions for complex growth rates are of primary 
interest. In the zero order we obtain the trivial solu- 
tion A^°^ = and in higher orders we find the nonzero 
corrections to A^"', corresponding (with accuracy up to 
the small terms of higher order) to two different types of 
deformation of the dust cloud borders Ci = ±C2: 

A^i) = ±ikdGce-^'"^, A(2) = ipfcdGcS (1 + 2kd) e'^'"^. 

(51) 

One solution corresponds to deformation of the cloud 
borders in the "phase," when (^i = (^2- This solution 
is stable since the real part of the growth rate is nega- 
tive. In the other solution, allowing for deformation of 
the cloud borders in the "anti-phase," when Ci = — C2, 
the real part of the growth rate is positive, and hence 
this solution is unstable. In both solutions the imaginary 
part of the growth rate is nonzero suggesting that the 
behavior of the system is oscillatory. 

Note, that the obtained result is true for any nontrivial 
values the parameters Gc and S. Thus, we conclude that 
the dust cloud of a small, but finite width is unstable 
with respect to short wavelength perturbations. 



3. Dust cloud of an arbitrary width 

Let us proceed to the discussion of general numerical 
results, obtained for a cloud of arbitrary width d, and 
consider the case S = 1. 

The neutral curves for a fixed value of the concentra- 
tion parameter Gc = 100 are presented in Fig. [7| The 
instability regions are above the curve 1 (between the 
curve 1 and the line k = 0) and under the curve 2 (be- 
tween the curves 2 and 3). Curve 3 is a straight line at 
which Gr — 9\/3S. This line defines the area of existence 
of a dust cloud and corresponds to a dust cloud of zero 
width. The cloud exists in the domain that is above this 
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FIG. 7: Neutral curves (line 1 and 2), plotted for Gc = 100, 
S = 1. Regions of unstable behavior are above the line 1 and 
under the line 2. Line Gr = 9\/3S (line 3) corresponds to a 
cloud of the vanishing width. 



line. As can be seen, there is a narrow range of values of 
Gr, where the system is stable. The results (|51f) . valid for 
a thin cloud in the short wavelength limit, are in agree- 
ment with those obtained from numerical solution to the 
general boundary value problem H3fci|l - (|42|l . However, it is 
obvious from Fig. [71 that the short wavelength instability 
is not "the most dangerous" : although in some range of 
governing parameters the flow is stable with respect to 
short waves, but yet for all possible values of the cloud 
width d there exist unstable two-dimensional perturba- 
tions in the form of transversal rolls with finite values of 
k. It is interesting to note, that the steady state in the 
case of a thin cloud is unstable with respect to perturba- 
tions with any wave number k larger than some critical 
value. 

The global minimum of the curve 1 and the global 
maximum of the curve 2 over fc > in Fig. are char- 
acterized by the critical values of the Grashof number: 
Grmi„, Grmaa;, which are reached at the values of the 
wave number kmim k^ax respectively. Let us discuss 
the dependence of these critical parameters on the con- 
centration parameter Gc. Figure |S1 gives the stability 
diagram on the plane (Gc, Gr). The stability region is 
under the line 1 and above the line 2. Note, that the 
curves 1 and 2 refer to the upper and the lower branches 
of the dependence c?(Gc), respectively (see Fig. 0J. This 
becomes important for distinguishing between stability 
regions for the upper and lower branches in the hystere- 
sis zone (Fig.lSJ, bounded by lines 4 and 5: the stability 
region of the upper branch is defined by the curve 1, for 
the lower branch it is defined by the curve 2. As before, 
the line 3 is a line determining the region of a cloud of 
vanishing width, at which Gr = OVSS. 

It should be noticed that the section of the straight line 
Gc — 0, at 9-\/3 < Gr < 496.3 also refers to a stability 
region. A piece of this line, corresponding to the values 
< Gr < 9^/3 is stable as well, however, the basic steady 
state in this case is completely free of particles, i.e. a dust 
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FIG. 8: Stability curves for the upper (line 1) and the lower 
(line 2) branches of the dependence d (Gc) on a plane (Gc, 
Gr) at S = 1. The stability regions are under the line 1 and 
above the line 2. Line 3 is the line Gr — 9\/3S; lines 4 and 5 
represent bounds of the hysteresis zone (see Fig-.^J. 



FIG. 10: Stability curves for the upper (lines 1 and 6) and 
the lower (line 2) branches of d(Gc) on a plane (Gc, Gr). 
The stability regions are under the line 1 and 6 and above the 
line 2. Line 3 is the line Gr = 9\/3S; lines 4 and 5 represent 
bounds of the hysteresis zone. 



cloud is absent. 

The variation of kmin, kmax with the concentration 
parameter Gc is shown in Fig. |3 The imaginary parts 
of growth rates corresponding to these solutions are 
nonzero, which means that this "concentration" mode 
is oscillatory. The hydrodynamic perturbations are not 
the most dangerous for our particular value of sedimenta- 
tion parameter S, their stability bound lies much higher 
than the lines 1 and 2. Of special note is the dependence 
^mm (Gc) at relatively small values of the concentration 
parameter Gc. It is seen, that for not very large values of 
Gc the wave number kmin is rather small, i.e. the system 
demonstrates a nearly long wavelength behavior. 

With the increase of the sedimentation parameter S the 
region of hysteresis wedge shifts to the region of larger 
values of Gr and Gc. Fig. EH presents the stability dia- 
gram for the case 8 = 5. The stability curves 1 and 2 do 
not undergo any qualitative changes. However, in con- 
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FIG. 9: Critical wave numbers kmin, kmax as functions of Gc 
at S = 1. 



trast to the case 8 = 1, the concentration mode is in com- 
petition with the hydrodynamic mode: there appears the 
range of values of Gc, at which the hydrodynamic per- 
turbations (line 6) become the most dangerous for the 
upper branch of the dependence d{Gc). With increase of 
S this range grows, reducing the total stability region. In 
the hysteresis zone line 6 sticks to the line 5: the upper 
branch becomes unstable here. In comparison with the 
case Gc = 0, the hydrodynamic mode is no longer mono- 
tonic. Indeed, the monotonic character of hydrodynamic 
perturbations is specified^ by the symmetry of the ve- 
locity profile uq. In our case, as it follows from H25|) . the 
symmetry of uq is broken at any value Gc ^ 0. 



4- Arbitrary three-dimensional perturbations 

It was mentioned before, that after linearization 
Eqs. H31() - (|33|l . boundary conditions (|34l) and the con- 
ditions at the borders of the dust cloud H29|) admit 
a transformation, which is analogous to the Squire 
transformation.'^ Under such a transformation the 
Prandtl number Pr does not change, whereas the param- 
eters Gr, Gc, and S transform according to the following 
rule: 



Gr 



Gr', Gc = -^Gc', S 



S', 



(52) 



where Gr, Gc, S are the parameters of the two- 
dimensional problem and Gr', Gc', S' are the correspond- 
ing parameters of the full three-dimensional problem; ky 
and kz are the wave numbers along the axes y and z, 
respectively; k-^ — ky + k^. Note, that at any given value 
kz/k the parameters of the two-dimensional problem Gr, 
Gc, S are always less than those of the three-dimensional 
problem. However, this fact does not imply that two- 
dimensional perturbations are the most dangerous. 
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Let us introduce the parameter, describing deviation 
of arbitrary three-dimensional perturbations from the 
two-dimensional perturbations in the form of transverse 
rolls H35|l . namely the angle between the wave vector 
k = (0, ky, kz), which generally lies in the plane y-z, and 
the z-axis. If we denote the cosine of this angle by a, 
then it follows from H52|) 

Gr' = — , Gc'^— , S' = -. 
a a a 

The greatest distinction of arbitrary three-dimensional 
perturbations from transversal rolls H35|l corresponds to 
the limiting case of longitudinal rolls, when a = (k^ = 
0, fc ^ 0). The perturbations of this kind are often called 
the "helical" perturbations, since the trajectory of a fluid 
element in such a flow looks like a helix. The motion 
of a fluid particle can be considered as a sum of two 
different forms of motion: first, the particle moves circle- 
wise inside the longitudinal roll, second, it is carried along 
the roll with the velocity uq + u. 

Thus, at given a the diagram of stability with respect 
to arbitrary three-dimensional perturbations can be ob- 
tained from that of two-dimensional perturbations l|35(l 
by rescaling Gr and Gc by a factor of l/a. It is important 
to note, that the parameter S is also transformed. 

The linear stability analysis, performed with respect to 
two-dimensional perturbations (|35|l for S from 1 up to 5, 
allows us to investigate the influence of three-dimensional 
perturbations with a in the range from 0.2 up to 1 for 
S = 5. It turns out, that in this range the most dangerous 
are the perturbations in the form of transversal rolls. 

Further, it can be explicitly shown, that the limit- 
ing case of the longitudinal rolls (helical perturbations), 
when a = (fcz =0, fc 7^ 0), docs not cause instability. 
Indeed, as it follows from Eqs. H31() - (|34() and linearized 
form of conditions H29|l . H3U|I . this particular case is de- 
scribed by the following boundary value problem (the 
equations are the same for each of three zones) : 

\u = —q + u" — kyii, (53) 

Xi) = ikyq + v" - kyi), (54) 

u' — ikyV = 0, (55) 

Xih + u'„u = w" - kyfv + Grd, (56) 

A^ + ^ = _L(^"_fc2^), (57) 

x = ±l: u = 0, w = 0, w = 0, 7? = 0, (58) 
X = a;i,2 : [q] = 0, [u] = 0, [v] = 0, [d] = 0, 
m = 0, N-0, K] = -GcCi,2, 
ACi^2 = u- (59) 

The problem for the fields ii and v splits off and can 
be treated independently. Excluding the pressure q from 
Eqs. (|53|I - H55|) and taking into account the continuity of 



ii and v at the points xi and X2, we obtain the following 
boundary value problem: 

\{v"~klv) = v'^ -2klv" + klv, (60) 
X = ±1 : V ^ Q, w' = 0. (61) 

Let us multiply Eq. 1)60(1 by the complex conjugate i* and 
integrate by parts across the layer. Taking into account 
(|61|l . after elementary calculations we obtain 




where (...) = J^^ . . . dx. Since the problem H60|l - H61|) 

describes perturbations in a quiescent viscous uniform 
fluid, it cannot give birth to instability: the growth rate 
is proved to be real and negatively defined; the perturba- 
tions monotonically decrease with time, and, therefore, 
one should look for a mode with u = 0. Hence, with 
account of incompressibility relation l|55|) , the boundary 
conditions H58() for u and continuity of this field at the 
points xi and X2 we obtain u — 0. 

Since u = 0, equation (|57|) corresponds to diffusion of 
in motionless fluid, and therefore, from analogous consid- 
eration we conclude that perturbations decrease. Hence, 
the mode with ■)? = is to be found. As it follows from 
the condition H59() for (i and (2 , the solution with A 7^ 
is possible only if = and C2 = 0- Consequently, w 
and its derivative are continuous at the points xi, X2, and 
the problem for w coincides with the problem for "d with 
the rescaled time tPr~^. Thus, for any values of govern- 
ing parameters the basic state is proved to be stable with 
respect to perturbations in the form of longitudinal rolls. 

V. CONCLUSIONS 

The interaction of the vortex buoyancy convective flow 
laden with particles of dust has been investigated in the 
case when the volume concentration of the solid phase is 
small. Under these conditions the particles are partially 
carried away by the flow, but due to sedimentation their 
velocity differs from the velocity of a fluid. At sufficiently 
intensive thermal convection some portion of the particles 
can be captured by the flow, which eventually results in 
formation of a cloud of dust. 

If mass concentration of the particles is rather small, 
the particles do not influence the flow and the study of 
the formation of a dust cloud is reduced to a kinematic 
problem. The size of the dust cloud captured by the con- 
vective vortex monotonically increases with the growth of 
the thermal Grashof number from some threshold value 
Gr*. At the values Gr < Gr*, the formation of a steady 
cloud is impossible: sooner or later all the particles set- 
tle on the bottom of a cavity. This critical value of the 
Grashof number is determined by the cavity shape and 
thermal boundary conditions and is also proportional to 
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the sedimentation parameter S. The calculations for the 
case of an infinite vertical layer heated from the sidewalk 
give Gr, = 9^38 . 

With the increase of mass concentration of the parti- 
cles their influence on the flow becomes signiflcant. This 
happens when relative variations of mass concentration 
are comparable with the Boussinesq nonisothermality pa- 
rameter (36. Generally, the growth of mass concentration 
results in the decrease of flow intensity. Energy of the 
flow is partly consumed on the particle motion against 
the gravity and on the motion of particles together with 
the descending flow. This energy is not returned to the 
flow completely because of the dissipative loss. The cal- 
culations show that the interaction of the particles with 
the convective flow can lead to the hysteresis phenomena. 
Different flow patterns can arise at the same values of the 
governing parameters: the modes, in which a compara- 
tively small dust cloud strongly suppresses the ascending 
flow and the modes developing at comparatively large 
cloud when the density inhomogeneities are low and the 
suppression effect is insigniflcant. 

The linear stability of the basic steady state in the form 
of a dust cloud in an infinite layer, heated from sidewalls, 
has been investigated. It is shown, that in a relatively 
narrow range of the values of governing parameters this 
state is stable. Based on the results of the stability anal- 
ysis, we can infer that in the typical case of a closed 



cavity the steady state with a relatively large dust cloud 
is most likely to have a larger range of stability. In the 
case of an infinite layer such a steady state is broken by 
long wavelength perturbations. Obviously, such pertur- 
bations cannot arise in a closed container of a finite size. 

The performed investigation demonstrates that the de- 
scribed effects of two-way interaction of fluid and parti- 
cles can exist in natural environment or be realized ex- 
perimentally. The values of the sedimentation parameter 
S of order 1, used in the linear stability analysis, corre- 
spond, for example, to the case of rather fine particles 
of dust with Tp ~ 10~^ cm, suspended under gravity 
in a gaseous medium filling a laboratory-scale container 
i ~ 1 cm {rp/L ~ 10"^, 5 ~ 10^ Ga ~ 10^). The 
developed theory is quite general and can be applied to 
describe similar effects not only in dusty media, but also 
in aerosols, liquids laden with small solid particles, and 
biological species in aqueous media. 
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